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The magnetization process, the susceptibility and the specific heat of the spin-1/2 AF-AF-F and 
F-F-AF trimerized quantum Heisenberg chains have been investigated by means of the transfer 
matrix renormalization group (TMRG) technique as well as the modified spin-wave (MSW) theory. 
A magnetization plateau at m = 1/6 for both trimerized chains is observed at low temperature. 
The susceptibility and the specific heat show various behaviors for different ferromagnetic and 
antiferromagnetic interactions and in different magnetic fields. The TMRG results of susceptibility 
and the specific heat can be nicely fitted by a linear superposition of double two-level systems, 
where two fitting equations are proposed. Three branch excitations, one gapless excitation and two 
gapful excitations, for both systems are found within the MSW theory. It is observed that the MSW 
theory captures the main characteristics of the thermodynamic behaviors at low temperatures. The 
TMRG results are also compared with the possible experimental data. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 



Physical properties of low-dimensional magnetic ma- 
terials have attracted considerable attention at the bor- 
der of condensed matter physics. In particular, the 
bond-alternating quantum magnets with both ferromag- 
netic (F) and antiferromagnetic (AF) couplings have 
been extensively investigated, because they have exhib- 
ited quite fascinating properties owing to the compe- 
tition of F and AF interactions. A simple AF bond- 
alternating chain (BAC) may be a F-AF alternating 
quantum Heisenberg chain, and the spin-1/2 copper- 
based compound (/P^4)CwC73pJ is a typical example, 
in which the Haldane-like behavior has been observed. 
This kind of alternating chains have been actively studied 
both theoretically and experimentally in the last decades. 
Another interesting antiferromagnetic BAC is the F-F- 
AF alternating Heisenberg chain with period n = 3. 
According to Oshikawa-Yamanaka- Affleck (OYA)|2|, the 
F-F-AF chain with spin S may exhibit a magnetiza- 
tion plateau at m = S/3, where m is the magnetiza- 
tion per site. In a typical example of spin-1/2 F-F-AF 
chain of 3CuCl 2 ■ 2dx with strong F coupling (Jf) and 
weak AF coupling ( Jaf)> however, no plateau in the low- 
temperature magnetic curve has been observed 3J . Such 
a breakdown of the magnetization plateau in small bond 
ratio of | Jaf/</f | El suggests that the OYA condition is 
only a necessary condition. 

There is another intriguing BAC with period 
n = 3, which consists of AF-AF-F alternations, 
and is a ferrimagnet. The examples of AF-AF-F 
chains are \Mn(L) 2 (N 3 ) 2 ] n with spin S = 5/20, 
[M(4,4'bimj)(N 3 ) 2 } n with M = Co (S = 3/2) and Ni 
(S = 1)0, and [Mn(N 3 ) 2 (bpee)] n with S = 5/20/. 



These compounds have strong F and weak AF cou- 
plings, and no magnetization plateaus have been seen 
experimentally. If the ratio \ Jaf/Jf\ are large enough, 
the plateaus could be expected during the magnetiz- 
ing process. We note that the plateau has been ob- 
served experimentally in a F-F-AF- AF BAC compound, 
Cu{3 - Cl P y) 2 (N 3 ) 2 11 M, which possesses a relatively 
large ratio | Jaf / Jf I \Ml LLJ] • 

Inspired by the exotic magnetic properties of the BAC 
compounds with period n = 3 observed experimentally, 
in this paper, we shall invoke the transfer matrix renor- 
malization group (TMRG) technique 01 to study the 
thermodynamic properties of the AF-AF-F chain and F- 
F-AF chain, respectively. The rest of this paper is or- 
ganized as follows. In Sec. II, we shall construct the 
model Hamiltonian for the trimerized J- J- J' Heisenberg 
spin chain. In Sec. Ill, we shall present our TMRG re- 
sults on the thermodynamic behaviors of the systems. 
The comparisons between the TMRG results and the re- 
sults of the modified spin-wave (MSW) theory will be 
discussed in Sec. IV. Finally, a brief summary will be 
given. 



II. MODEL 

Let us consider a trimerized S — 1/2 J- J- J' Heisen- 
berg quantum spin chain. The Hamiltonian of the system 
reads 

H = ^](JS3j_2 • S3j_i + JS3j_i • S 3 j + J'S 3 j ■ S3J + 1) 
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FIG. 1: Spin arrangements of the J-J-J' trimerized spin chain 
with 5* = 1/2, where | denotes spin up, j denotes spin down, J 
and J' are exchange couplings, (a) AF-AF-F chain, J = Jaf, 
J' = J F ; (b) F-F-AF chain, J = J F , J' = J AF . 

where J and J' are exchange integrals with J, J' > de- 
noting the antiferromagnetic coupling and J, J' < the 
ferromagnetic coupling, h is the external magnetic field, 
and we take g^iB = 1- The schematic spin arrangements 
of the system are shown in Figs. When J = J' = J F 
(or Jaf), the system becomes a uniform 5=1/2 Heisen- 
berg ferromagnetic (or antiferromagnetic) chain. 

When J ^ J', as indicated in Figs. Q if J = Jaf, 
J' = Jp, it is a ferrimagnetic spin chain (i.e. the config- 
uration is AF-AF-F); if J = Jp, J' = Jaf, it becomes 
an antiferromagnetic system (i.e. the configuration is F- 
F-AF). The competition of F and AF couplings in these 
two trimerized chains would give rise to plenty of inter- 
esting characteristics. The magnetic properties of this 
model with S=l/2, 1, 3/2 and 2 in a magnetic field at 
zero temperature have been investigated by use of the 
density matrix renormalization group method [l2T|. It is 
purpose of this paper to study the thermodynamic prop- 
erties of these trimerized quantum Heisenberg spin chains 
by means of the TMRG technique 01 . 

III. TMRG RESULTS 

The TMRG method was detailed in two nice 
reviews 14f, and we shall not repeat the technical skills 
here for concise. In our calculations, the number of kept 
optimal states is taken as to = 64 for the susceptibility 
X, and to = 80 for the specific heat C, where the width 
of the imaginary time slice is taken as e — 0.1. We have 
used different to and e to verify the accuracy of calcu- 
lations. At high temperature, the error caused by the 
Trotter-Suzuki decomposition is important and is of the 
order e 3 for a fixed to. At low temperature, the error 




h/J F T/J F 

FIG. 2: (Color online) For the AF-AF-F chain with 
Jaf/Jf = 1: (a) The magnetization per site m as a func- 
tion of magnetic field h at different temperatures; (b) The 
temperature dependence of m under different magnetic fields. 



resulting from the basis truncation becomes important, 
which drops exponentially with increasing to initially and 
reaches to a finite value. The Trotter-Suzuki error is less 
than 10 -3 , and the truncation error is smaller than 10~ 6 
in our calculations. The physical quantities presented 
below are calculated down to T — 0.025 (in units of Jp). 



A. MAGNETIZATION 

Fig. HJa) shows the magnetization per site to as a 
function of magnetic field h at different temperatures 
for the AF-AF-F chain with Jaf/Jf = 1- At zero 
temperature, a magnetization plateau at to — 1/6 was 
observed^. As illustrated in Fig. Qfa), the ground 
state of the 5=1/2 AF-AF-F chain in the absence of 
the magnetic field corresponds to to — 1/6. When the 
external field is applied, the ground state with to = 1/6 
still persists for h/Jp < 1.3, which is consistent with 
the OYA condition 0- For hj Jp > 1.3 the plateau state 
is destructed, and m(h) increases rapidly with increasing 
the field. When the field is increased to hj Jp = 1.6, m(h) 
becomes saturated, and the ground state becomes a fully 
polarized state. When the temperature is increased, the 
magnetization plateau retains at low temperature, and is 
gradually smeared out at high temperature owing to the 
thermal fluctuations. It is noting that the magnetization 
m(h) as a function of h at finite temperature starts from 
zero, that is nothing but the result of Mermin- Wagner 
theorem. 

Fig. Efb) gives the temperature dependence of the 
magnetization m(T) under different magnetic fields for 
the AF-AF-F chain with Jaf/Jf = 1- It can be seen 
that to(T) behaves distinctly for different ranges of the 
external field. At ft, < 0.6, to(T) first decreases rapidly, 
and then declines slowly with increasing temperature; at 
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h/J F T/J F 

FIG. 3: (Color online) For the F-F-AF chain with Jaf/Jf = 
1: (a) The magnetization per site m as a function of the 
magnetic field at different temperatures; (b) The temperature 
dependence of m under different fields. 



0.6 < ft < 1.4, m(T) first increases, and then declines 
slowly with increasing temperature; at ft > 1.4, m(T) de- 
clines with increasing temperature. At ft = 0.6 and 1.4, 
m(T) appears to remain constant with zero curvature at 
low temperature, and declines slowly with temperature. 
It is found that just below and above h = 0.6 and 1.4, 
the curvature of m(T) changes abruptly with opposite 
signs, implying that ft = 0.6 and 1.4 could be viewed as 
the crossover fields. This observation is also consistent 
with Fig. |2Ia) where there are two crossing points in the 
curves for different temperatures at ft — 0.6 and 1.4. 

Fig. Efa) shows the field dependence of the magneti- 
zation per site m(h) at different temperatures for the F- 
F-AF chain with Jaf/Jf = 1. A magnetization plateau 
at m = 1/6 is also observed at zero temperature. As 
indicated in Fig. Hfb), the ground state of the S = 1/2 
AF-AF-F chain in the absence of a magnetic field corre- 
sponds to m — 0. With increasing the magnetic field, m 
first increases till to m = 1/6 at h/Jp ~ 0.2, then goes 
into the plateau state with m — 1/6 for h/Jp < 0.7, be- 
yond which the plateau state is destroyed, and m(h) then 
increases rapidly to a saturation at h/Jp — 0.8 where 
the system becomes fully polarized. The occurrence of 
the plateau state with m = 1/6 is also consistent with 
the OYA condition|2|. It is seen that at low temperature 
the magnetization plateau still persists, though the width 
of plateau is getting narrow. At high temperature, the 
plateau is gradually smeared out by thermal fluctuations, 
and m(h) increases nonmonotonically till the saturation 
with increasing the magnetic field. 

The temperature dependence of the magnetization per 
site m(T) under different fields for the F-F-AF chain with 
Jaf/Jf = 1 is presented in Fig. |3Jb). One may see that 
m(T) shows different behaviors at low temperature in 
different ranges of the magnetic field. At ft, < 0.2, m(T) 



increases first and then declines smoothly with increasing 
temperature; at 0.2 < ft < 0.4, m(T) decreases dramat- 
ically at low temperature and then slowly with increas- 
ing temperature; at 0.4 < ft < 0.7, m(T) increases at 
low temperature and then declines slowly with increasing 
temperature; at ft > 0.7, m(T) declines remarkably with 
increasing temperature. At ft, = 0.4 and 0.7, m(T) re- 
mains almost unchanged at lower temperature, and then 
declines slowly with increasing T. Again, the curvature 
of m(T) changes abruptly with opposite signs just be- 
low and above ft = 0.4 and 0.7. The results presented 
in Fig. OUb) suggest that at low temperature the system 
may enter into different states in different ranges of the 
magnetic field, namely, for ft < 0.2, the system is in an 
ordering state; for 0.2 < ft < 0.4 the system could go 
into the spin-canting state; for 0.4 < ft < 0.7 the sys- 
tem enters into the plateau state where the excitations 
are gapful; and for ft > 0.7 the system goes into another 
ordering state before saturation. 

It is interesting to note that the magnetization as a 
function of magnetic field and temperature exhibits dif- 
ferent behaviors for the AF-AF-F chain and the F-F-AF 
chain, for the former is a ferrimagnet, while the latter 
is an antiferromagnet. However, both systems shows a 
magnetization plateau with m = 1/6 at low temperature, 
being in agreement with the OYA condition. 



B. SUSCEPTIBILITY 

Let us now present the TMRG results of the suscepti- 
bility (x) as a function of temperature (T) for the AF- 
AF-F chain. In the absence of the magnetic field, the sus- 
ceptibility diverges when the temperature tends to zero, 
and with increasing temperature, x decreases rapidly, as 
shown in Fig. 0£a). For different Jaf/Jf, x(T) ex- 
hibits similar behaviors at low T, but distinct behaviors 
at moderate T, as manifested in the inset of Fig. Ufa), 
where \T versus T is plotted for various Jaf/Jf, and 
a low-T divergence and a round minimum are observed 
for Jaf/Jf < 1, indicating that when the F interaction 
Jf predominates, the system shows a ferrimagnetic be- 
havior; when the AF interaction J af is predominant, the 
system falls into a trimerized bond state. In the presence 
of the magnetic field, the susceptibility x shows quite 
different behaviors, as presented in Figs. Efb) and (c) at 
Jaf/Jf = F Unlike the case of ft = 0, x under nonzero 
ft exhibits a peak at low temperature except ft = 0.6 and 
1.4 where with increasing T, x first increases and then 
drops slowly for the former, and goes to divergence at 
low temperature for the latter. With the increase of ft, 
the peak of x appears to become rounded, and the peak 
height of x gets suppressed for ft < 0.6, then increased 
for 0.6 < ft < 1.4, and then becomes suppressed again for 
ft > 1.4. These observations show again that the system 
can enter into different states within different ranges of 
the magnetic field. 

For the F-F-AF chain, the temperature dependence 
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FIG. 4: (Color online) The temperature dependence of the 
susceptibility \. For the AF-AF-F chain (a) at h = for 
different Jaf/Jf, (b) and (c) at Jaf/Jf = 1 for different 
external fields, where the inset of (a) depicts \T against T 
for different Jaf/Jf- For the F-F-AF chain (d) at h = 
for different Jaf/Jf, (e) and (f) at Jaf/Jf = 1 for different 
external fields. 
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FIG. 5: (Color online) (a) The experimental data of \T of 
the spin-5/2 AF-AF-F chain [Mm{N 3 ) 6 {bpe)3\ (from Ref.Q); 
Inset presents the TMRG result of an AF-AF-F chain with 
S = 1/2. (b) A comparison of the susceptibility for SCuCh ■ 
2dx with the TMRG result of a F-F-AF chain with S = 1/2, 
where the experimental data are taken from Ref. Q. 



vanishes at h/Jp = 0.3, while the second round peak 
becomes slightly higher, as shown in Fig. E{e). When 
h/Jp > 0.5, the second peak of x moves to the low- 
temperature side, and at h/Jp = 0.7, the two peaks 
merge into a single; when h is increased further, the peak 
of x is suppressed again. The behaviors of x(T) show that 
the F-F-AF chain can enter into different states at low 
temperatures under different magnetic fields, consistent 
with the results presented in Fig. E^b). 



of the susceptibility x f° r different Jaf/Jf and h is 
shown in Figs. Efd)-(f), respectively. In the absence 
of the magnetic field, the susceptibility reveals peaks 
at low temperature for different Jaf/Jf, and x = 
at T = 0. When Jaf/Jf is either smaller or larger, 
e.g., Jaf/Jf = 0.25, 2, 4, x shows one peak; while 
Jaf/Jf = 0.5 and 1, x shows a double-peak structure, 
as manifested in Fig. EJd). The calculated results indi- 
cate that the system shows an AF behavior. It seems 
that the double-peak structure comes from the two dif- 
ferent kinds of excitations induced by the competition 
between F and AF interactions. In the presence of the 
magnetic field, the temperature dependence of the sus- 
ceptibility x at Jaf/Jf = 1 for different external fields 
are presented in Figs. Hfe) and (f). It is found that x(T) 
behaves variously in different ranges of the magnetic field. 
When h is small, x(T) shows a double-peak structure at 
low temperature with the first peak sharp and the sec- 
ond peak round; with increasing h, the height of the first 
peak grows dramatically, while that of the second round 
peak leaves almost unchanged; when h/Jp > 0.2, the 
first peak is remarkably suppressed at h/Jp = 0.25, and 



To verify our TMRG results, we have included 
the experimental data of an AF-AF-F chain com- 
pound \Mn^{N^)Q{bpe)z} and a F-F-AF chain compound 
2>CuCh ■ 2dx in Figs. C3 for a comparison. Fig. 

shows the experimental data of xT versus T for 
[Mn3(N3)e(bpe)s], whose magnetic structure is an AF- 
AF-F chain with spin-5/2 |7j. Owing to our comput- 
ing capacity, we cannot at present calculate directly the 
magnetic properties of the spin-5/2 AF-AF-F chain by 
using the TMRG method. As a qualitative understand- 
ing, however, in the inset of Fig. ISTa). we have pre- 
sented our TMRG result of a spin- 1/2 AF-AF-F chain 
with h/Jp — 0.003 and Jaf/Jf — 0.2. It can be seen 
that the experimental and calculated curves share quali- 
tatively similar characteristics. Fig. \Ej[b) gives a com- 
parison of the experimental susceptibility of the com- 
pound ?>CuCli ■ 2dx, which is a F-F-AF chain with spin- 
1/2 |, with the TMRG result. It is found that our 
TMRG result fits well with the experimental data, with 
the coupling parameters Jaf/Jf — 0.2, Jp — 100 K, 
and h/Jp = 0.005. Our calculated result is in agreement 
with the previous theoretical result |3j. 
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FIG. 6: (Color online) The temperature dependence of the 
specific heat C. For the AF-AF-F chain (a) at ft = for 
different Jaf/Jf, (b) and (c) at Jaf/Jf = 1 for different 
external fields. For the F-F-AF chain (d) at h — for different 
Jaf/Jf, (e) and (f) at Jaf/Jf = 1 for different external 
fields. 



C. SPECIFIC HEAT 

In this subsection, we shall discuss the specific heat 
of the spin- 1/2 AF-AF-F and F-F-AF chains, as shown 
in Figs. El for different Jaf/Jf and h. For the AF- 
AF-F chain, in absence of the magnetic field, the spe- 
cific heat C as a function of temperature shows clearly a 
double-peak structure, as indicated in Fig. Et a )- With 
increasing Jaf/Jf, the low-temperature peak of C leaves 
almost unchanged, while the second peak moves to the 
high-temperature side, and gets rounded. This feature of 
the specific heat shows that the system might exist two 
different kinds of excitations due to the competition of 
F and AF interactions, as will be discussed in next sec- 
tion. It seems that the first peak corresponds to the F 
excitations, while the second peak corresponds to the AF 
excitations. 

In presence of the magnetic field, the temperature de- 
pendence of the specific heat C for the AF-AF-F chain 
with Jaf/Jf = 1 is shown in Figs. EJb) and (c) under 
different fields. When h is small, the double-peak struc- 
ture of C is observed; with increasing h, the first peak 
tends to vanish, while the second peak remains almost 



unchanged, as shown in Fig. EJb). When 0.2 < h/Jp < 
1.1, the two peaks become a single. When h/Jp > 1.1, 
as presented in Fig. Etc), the first peak recovers anew, 
and with the increase of h, the first peak moves to the 
lower temperature side till h/Jp = 1.4, then shifts to the 
high temperature side for 1.8 > h/Jp > 1.4, and tends 
to vanish for h/Jp > 1.8, while the second peak shifts 
slightly to the higher temperature side with increasing 
h. The results show that the system is in a few different 
states under different magnetic fields. 

For the F-F-AF chain, the temperature dependence of 
the specific heat C is shown in Figs. Etd)-(f) for different 
Jaf/Jf and h. In absence of the magnetic field, C also 
shows a double-peak structure; with increasing Jaf/Jf, 
the height of the first peak of C is suppressed but its 
position leaves unchanged, while the second peak tends 
to shift to the high-temperature side, as illustrated in Fig. 
EJd). It is seen that the qualitative behavior is similar 
to that of the AF-AF-F chain presented in Fig. Eta). 
The results suggest that the F-F-AF chain has also two 
kinds of excitations due to the competition of F and AF 
interactions. 

Figs. Et e ) and (f) gives the temperature dependence of 
the specific heat C for the F-F-AF chain with Jaf / Jf — 
1 under different magnetic fields. Similar to the AF-AF- 
F chain, with increasing h, the first peak of C is gradually 
suppressed, while the second peak retains nearly intact 
till h/Jp = 0.2. When 0.2 < h/J F < 0.5, the first peak 
is completely suppressed; when 0.5 < h/Jp < 0.7, with 
the increase of h, the first peak recovers again with the 
peak height decreased and the peak position shifted to 
the lower temperature side, while the second peak moves 
slightly to the high temperature side with the peak height 
a little bit enhanced; when h/Jp > 0.7, the situation 
changes, namely, with increasing h, the first peak of C is 
enhanced and moves to the high temperature side, and 
the second peak is also promoted. The present obser- 
vation displays that the F-F-AF chain falls into differ- 
ent thermodynamic states under different external fields, 
and has some behaviors similar to the AF-AF-F chain, 
although both chains have quite different ground states, 
as the latter is a ferrimagnet, while the former is an an- 
tiferromagnet. 



D. FITTING TO A SUPERPOSITION OF 
DOUBLE TWO-LEVEL SYSTEMS 

The double-peak structure of the temperature depen- 
dence of the specific heat can be viewed as a minimum at 
low temperature and a Schottky-like maximum at high 
temperature. The Schottky-like anomaly can be fitted 
with a so-called two-level system, where the thermal pop- 
ulation of the levels is governed by the Maxwell-Bolzman 
statistics, n, ; = e~ Ei l kBT /Z (i = 1,2) with the distribu- 
tion function Z = e - £ i/ fc s T + e -~^/k B T ^ and the total 

energy E(T) — n\E\ + n 2 e 2 - The specific heat can be 
obtained by C = {dE/dT). Define ei = 0, e 2 = A. The 
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Schottky-like specific heat can be written as 



C a (A,T) = k B ( 



= A/fc B T 



A 



k B T> (l + e A/fc B T)2 



(2) 



Similarly, it is found that the magnetic susceptibility can 
be fitted by [3, 



Xs (A,T) = k B ( — -) 



„A/k B T 



k B T> (l + e A / fc « T ) 2 ' 



(3) 



where A is the excitation gap. However, our TMRG 
results cannot be well fitted by these above two formulae. 
By noting that the present systems such as the AF-AF-F 
and F-F-AF chains are quantum trimerized systems with 
a period n — 3, and the double-peak structure may be 
regarded as the consequence of the competition between 
F and AF excitations, the thermodynamic properties of 
these systems could be mimicked by a superposition of 
double two-level systems. As a result, the susceptibility 
for the AF-AF-F and F-F-AF chains may be expected to 
fit by 



X (T) = A lXs (A u T) + A 2Xs (A 2) T), 



(4) 



and the double-peak structure of the specific heat can be 
fitted by the following form: 



C{T) = AiC s (Ai,T) + A 2 C S (A 2 ,T), 



(5) 



where A\, A 2 are fitting parameters, and A\, A 2 are the 
excitation gaps for the double two-level systems. In these 
above equations, Xs and C s are defined by Eqs. © and 
©, respectively. 

For the AF-AF-F chain, in absence of the magnetic 
field, we have found that x(T) cannot be fitted by Eq. 
(@J, but can be well fitted by the equation x(T) — 
a{T/J F )- 2 - b{T/J F y 3 / 2 + ciT/Jp)- 1 with a, 6, and 
c the fitting parameters, as shown in Fig. Ufa) for 
Jaf/Jf — 0.5 as an example. We note that the sus- 
ceptibility of a spin- 1/2 ferromagnetic system was ex- 
panded as xJ/L{gii B ) 2 = 0.04167i~ 2 - 0.145*- 1 - 5 + 
O.m- 1 + 0(t-°- 5 ) with t = —k B T I J (J < 0) at low 
temperatures [T^ . Thus, it is not surprising that the 
present system shows a similar behavior, as the AF-AF- 
F chain is a ferrimagnet. However, in presence of the 
magnetic field, the susceptibility x{T) of the AF-AF-F 
chain can be nicely fitted by Eq. as compared in 

Fig. Efb) for Jaf/Jf — 1 and h/Jp — 1.7 as an exam- 
ple. Figs. E[c) and (d) give the TMRG fitting results 
in terms of Eq. Q for the susceptibility x(T) of the 
F-F-AF chain in absence and presence of the magnetic 
field, respectively. One may see that apart from the case 
where x(^) 1S divergent at lower temperature and cannot 
be fitted by Eq. (0J, the susceptibility as a function of 
temperature for both trimerized chains can be well fitted 
by a superposition of the double two- level systems. 

Figs. |H1 present the TMRG results of the specific heat 
as a function of temperature fitting to a superposition 
of double two-level systems characterized by Eq. (J5J for 
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FIG. 7: (Color online) The TMRG results of the susceptibility 
as a function of temperature are fitted. For the AF-AF-F 
chain (a) by a(T/ J F y 2 -b(T/ J F y 3/2 +c(T/ Jp)" 1 for h = 0, 
(b) by Eq. J1J for h nonzero; for the F-F-AF chain by Eq. 
|@0 for (c) h — 0, and (d) h nonzero. 
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FIG. 8: (Color online) The TMRG results of the specific heat 
as a function of temperature are fitted by Eq. (|KJ. For the 
AF-AF-F chain with (a) Jaf/Jf = 0.5 and h = 0, and (b) 
Jaf/Jf = 1.7 and h/J F = 1.7; for the F-F-AF chain with 
Jaf/Jf = 1 and (c) h = and (d) h/J F = 0.85. 
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the AF-AF-F and F-F-AF chains. It can be found that 
the specific heat with double-peak structure can be well 
fitted by Eq. J^J at low temperature, with only a slightly 
quantitative deviation at high temperature, showing that 
the main features of the specific heat of both trimerized 
chains can be reproduced by a superposition of double 
two-level systems. The fitting results could give the two 
excitation gaps. The reasons for such nice fittings come 
from the fact that these two trimerized chains have three 
branches of excitations (see next section) , one gapless F 
excitations, and two gapful excitations, which could be 
equivalently treated by a superposition of double two- 
level systems with two gaps Ai and A2. We have also 
found that a direct fitting to a three-level system is not 
as good as to double two-level systems. 



IV. MODIFIED SPIN- WAVE THEORY 

The zero-field specific heat of a quantum ferrimagnet, 
which is believed as an intrinsic double-peak structure 
of topological origin and is different from the external 
field induced double-peak structure, was inv estigate d b y 
using the modified spin- wave (MSW) theory jl7t 1 18l Il9f . 
where it has been found that in the case of the AF-AF- 
F-F BAC, there are four distinct branches of spin- wave 
excitations, say, the lower two ferromagnetic bands con- 
struct the low-temperature bump, while an upper ferro- 
magnetic band and an antifcrromagnetic band contribute 
to a Schottky-type peak at mid temperatures. 

The field-induced double-peak structure of the specific 
heat in a quantum ferrimagnet has also been studied by 
using the linear spin- wave theorvpol|. In the case of spin- 
(1, i) ferrimagnet, a simple picture was given by noting 
that the Zeeman term introduces a gap to the F excita- 
tions that increases with increasing the field, and a gap 
to the AF excitations that decreases with increasing the 
field, and the specific heat really reflects the dual struc- 
ture of the antiferromagnetic and ferromagnetic excita- 
tions for all fields. 

The MSW theory can also be applied to the present 
5 = 1/2 AF-AF-F and F-F-AF chains in absence of the 
external field so that our TMRG results of the thermo- 
dynamic quantities may be used to compare with the 
corresponding MSW results. In the conventional spin- 
wave scheme, the spin deviations in each sublattice di- 
verge in the one-dimensional (ID) antiferromagnets, but 
the quantum as well as thermal divergence of the number 
of bosons can be overcome in the Takahashi scheme [l7j 
that will be applied to the present antiferromagnetic F- 
F-AF chain. The AF-AF-F chain is a ferrimagnet, whose 
magnetization should be nonzero in the ground state but 
zero at finite temperature, leading to that we can apply 
the Yamamoto scheme pit |22| . where the Lagrange mul- 
tiplier was introduced directly in the free energy, to our 
present ferrimagnetic AF-AF-F chain. The detail deriva- 
tions of the MSW formalism are collected in Appendix 
A, where the linear modified spin-wave (LMSW) theory, 
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FIG. 9: (Color online) The excitation spectra of the AF-AF- 
F chain calculated by the LMSW and PIMSW methods with 
(a) Jaf/Jf = 0.25 and (b) Jaf/Jf = 4. A comparison 
of the TMRG results with the LMSW and PIMSW results 
of the susceptibility x(T) for the AF-AF-F chain with (c) 
Jaf/Jf = 0.25 and (d) Jaf/Jf = 4. The comparison of 
TMRG with LMSW and PIMSW results of the specific heat 
C(T) for the AF-AF-F chain with (c) Jaf/Jf = 0.25 and (f) 
Jaf/Jf = 4. 



which is up to the order of 0(S 1 ), and the perturbational 
interacting modified spin-wave (PIMSW) theory, which 
is up to the order of O(S ), are included. Our results 
show that the MSW results may describe well the low- 
temperature behavior of the system, but they are not 
good in agreement with the TMRG results at high tem- 
perature. 

Figs. OJa) and (b) present the zero-field excitation 
spectra which were calculated by the LMSW and PIMSW 
methods within the Yamamoto scheme for the AF-AF- 
F chain. It is seen that this trimerized system has 
three branches of excitation spectra, one gapless excita- 
tion spectrum and two gapful excitation spectra. When 
Jaf/Jf is small, the lower F gapless spectrum and the 
upper AF gapful spectrum look coincidence, while the 
mid F branches are slightly shifted; when Jaf/Jf be- 
comes large, the lower F branch is far separated from 
the mid and upper excitation spectra for the PIMSW 
results, while for the LMSW results, the excitation spec- 
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tra are more dispersive, and the mid branch has two 
crossing points with the upper branch, showing that the 
LMSW and PIMSW methods give different results when 
the AF interaction is more dominant. Figs. Et c ) and (d) 
show the comparison between the TMRG and LMSW 
and PIMSW results of the temperature dependence of 
the zero-field susceptibility for the AF-AF-F chain with 
different Jaf/Jf- When Jaf/Jf is large, the LMSW 
and PIMSW results of x(T) are in agreement with the 
TMRG result, as shown in Figs. Eld); when Jaf/Jf is 
small, say, the F interaction takes predominant, the low- 
temperature behavior of x(T) calculated by the LMSW 
and PIMSW theories coincides with the TMRG result, 
while the high-temperature behavior shows a somewhat 
difference, as shown Figs. Etc). Figs. Et e ) and (f) give 
the comparison of the LMSW and PIMSW results with 
the TMRG result of the zero-field specific heat C(T) for 
the AF-AF-F chain with different Jaf / Jf- It can be seen 
that the LMSW and PIMSW theories can be applied to 
describe the specific heat of the AF-AF-F chain at lower 
temperature, where the PIMSW result can fairly recover 
the first peak position of C(T) at lower temperature, but 
they cannot describe well the high-temperature behavior 
where both MSW results look only qualitatively similar 
to the TMRG result. However, it can be stated that the 
double-peak structure of the specific heat is intimately 
related to the three branch excitations of the trimerized 
system, as manifested by the MSW results. 

Figs. HUf a) and (b) present the excitation spectra cal- 
culated by the LMSW and PIMSW methods within the 
Takahashi scheme for the F-F-AF chain with Jaf/Jf = 
0.25 and 4 in absence of the magnetic field, respectively. 
There are also three branch excitations, i.e., one F gap- 
less excitation spectrum, and two gapful excitation spec- 
tra, which are responsible for the thermodynamic be- 
haviors of the system at low temperature. When the 
AF interaction is weak, the LMSW and PIMSW theo- 
ries give a similar result, as shown in Fig. HOf ah when 
the AF interaction becomes stronger, both methods pro- 
duce quite different results particularly for the upper ex- 
citation branches, say, the two upper excitation spectra 
calculated from the PIMSW method have two crossings, 
while those calculated by the LMSW method are sep- 
arated, as indicated in Fig. HUf bL showing that the 
LMSW and PIMSW methods should be cautiously ap- 
plied when the F interaction of the system is weaker. 
Figs. HUT c) and (d) show a comparison between the 
TMRG and MSW results of the temperature dependence 
of the zero-field susceptibility for the F-F-AF chain with 
Jaf/Jf — 0.25 and 4, respectively. One may see that 
when Jaf / Jf is small, the low-temperature peak of x(T) 
can be fairly recovered, although there are somewhat dif- 
ferences; when Jaf/Jf is large, say, the F interaction is 
weak, the position of low-temperature peak of x(T) can 
be clearly reproduced, although there are quantitative 
changes, indicating that the LMSW and PIMSW meth- 
ods could capture the main characteristics of the low- 
lying excitations of the present trimerized system at low 
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FIG. 10: (Color online) The excitation spectra of the F-F-AF 
chain calculated by the LMSW and PIMSW methods with 
(a) Jaf/Jf = 0.25 and (b) Jaf/Jf — 4. A comparison of 
the TMRG results with the LMSW and PIMSW results of the 
susceptibility x(T) f° r the F-F-AF chain with (c) Jaf/Jf = 
0.25 and (d) Jaf/Jf = 4. The comparison of TMRG with 
LMSW and PIMSW results of the specific heat C(T) for the 
F-F-AF chain with (e) Jaf/Jf = 0.25 and (f) Jaf/Jf = 4. 



temperature. Figs. HOf e) and (f) give the comparison be- 
tween the TMRG and MSW results of the zero-field spe- 
cific heat C(T) for the F-F-AF chain with Jaf /Jf = 0.25 
and 4, respectively. It is seen that the double-peak struc- 
ture of C(T) can also be obtained by the LSMSW and 
PIMSW methods, which can produce the main features 
of the specific heat of the system at low temperature, but 
have large deviations at high temperature. 



V. SUMMARY 

In this paper, we have numerically studied the thermo- 
dynamics of the S = 1/2 AF-AF-F and F-F-AF trimer- 
ized quantum Hciscnbcrg chains by means of the TMRG 
method and the MSW theory. The magnetization pro- 
cess, the low-lying excitations, the susceptibility and the 
specific heat of both systems have been explored. It is 
found that the TMRG results of the temperature depen- 
dence of the susceptibility as well as the specific heat can 
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be nicely fitted by a superposition of double two-level 
systems, and two fitting equations for the peak struc- 
tures are proposed. For the AF-AF-F chain with strong 
F interactions, x(^) can be well fitted by a polynomial. 

For the AF-AF-F chain which is a ferrimagnet, it is ob- 
served that the magnetic curve m(h) reveals some pecu- 
liar behaviors where a magnetization plateau at m = 1/6 
is observed at low temperature T, while the plateau is 
smeared out when T is increased, and the temperature 
dependence of the magnetization m(T) shows various be- 
haviors in different ranges of the magnetic field. The sus- 
ceptibility x(T) exhibits different behaviors in different 
fields: at h = 0, x(T) goes to divergent at lower T and 
decreases with increasing T for Jaf/Jf < 1; at /i ^ 0, 
x{T) shows a low-temperature peak for Jaf/Jf = 1 ex- 
cept h/Jp = 0.6 and 1.4, and decreases at high T, where 
the peak height and position of x(T) vary under different 
h. The behavior of \T against T for a spin-5/2 AF-AF- 
F chain compound [Mnz{N^)Q{bpe)-i\ looks qualitatively 
similar to the TMRG result of a spin- 1/2 AF-AF-F chain. 
The specific heat C(T) is found to show a double-peak 
structure, with the peak height and position varying for 
different Jaf/Jf as well as under different h. The results 
of m(T), x(T) and C(T) show that the system under in- 
terest could be in several different states under different 
magnetic fields. 

For the F-F-AF chain which is an antiferromagnet, it is 
found that the magnetization process shows an interest- 
ing behavior, and a magnetization plateau at m = 1/6 
is also seen at low temperature which is smeared out 
at high temperature. m{T) displays different behaviors 
in different regimes of the magnetic field, which is fur- 
ther confirmed by the TMRG results of the susceptibil- 
ity x(T) where a peak is observed at low temperature, 
but the peak height and position are found dependent 
on different ranges of h. The experimental data of x(T) 
for the spin-1/2 F-F-AF chain compound ZCuCl^ • 2dx is 
compared with the corresponding TMRG result, which 
is found quite agreement. The specific heat as a function 
of temperature shows a double-peak structure, whose 
shapes change with different Jaf / Jf as well as in dif- 
ferent magnetic fields. 

The MSW theory within the Yamamoto scheme as 
well as the Takahashi scheme is also applied to study 
the thermodynamics of the present trimerized quantum 
Heisenberg AF-AF-F chain and F-F-AF chain, respec- 
tively. Three branch excitations are found for both sys- 
tems, say, one gapless excitation spectrum and two gap- 
ful excitation spectra, that come from the competition 
between the F and AF interactions. A perfect fitting of 
x(T) and C(T) to a superposition of double two-level sys- 
tems could be reasonably understood within the frame- 
work of the MSW theory. It is three branch excitations 
with two gaps in the long wave-length limit that are re- 
sponsible for the peak structures of the susceptibility and 
the specific heat observed by the TMRG calculations. It 
is found that the MSW theory could to some extent de- 
scribe the low-temperature behavior of the systems, and 



gives only a qualitative description for x(T) and C(T) 
of AF-AF-F chain and x(T) of F-F-AF chain at high 
temperature, but has a large deviation for C(T) of the 
F-F-AF chain at high temperature. 

Finally, we would like to state that the present TMRG 
and MSW results could be helpful for understanding 
the magnetization process as well as the thermodynamic 
properties of the AF-AF-F and F-F-AF trimerized quan- 
tum spin chains. As there are few experimental data 
available for such systems at present, we expect that more 
experiments can be performed in near future to verify our 
theoretical findings. 

APPENDIX A: MSW THEORY FOR THE 
TRIMERIZED F-F-AF CHAIN 

The formalism of the MSW theory for the AF-AF-F 
chain can be found in Ref. 23]. In this Appendix, the 
MSW theory for the F-F-AF chain is formulated. The 
Hamiltonian of the ID spin-1/2 F-F-AF chain can be 
written as 

N 

H = V^(JpSl„ • S2 n + Jf^2u ■ S3n + JAF^3n ■ S4 n 

71=1 

J F S + J F s + JafS 6 „ • Si„ + i), 

where 6N is the length of the F-F-AF chain and Jf < 0, 
Jaf > 0. We start from the Holstein-Primakoff (HP) 
transformation 

Sj n = —S + a~^ n a,j n , (Al) 

where S = 1/2, i = 1, 2, 3 and j = 4, 5, 6. The LMSW 
method treats the HP transformation up to 0(S 1 ), while 
the PIMSW method treats is the HP transformation up 
to O(S ) where the interactions between spin waves are 
handled perturbatively. 

In the conventional spin-wave scheme, the spin devia- 
tions in each sublattice, (af k aik), diverge in the antiferro- 
magnetic ground state but stay finite in the ferrimagnetic 
ground state. In the AF system, the quantum as well as 
the thermal divergences of the number of bosons are all 
suppressed: 

n 

6 

= 6YS-^5>+a lfc >. (A2) 

k i=l 

As the F-F-AF is an antiferromagnet, the Takahashi 
scheme 17] will be invoked. In order to enforce the con- 
straint (|A2|) , a Lagrange multiplier fi is introduced to the 
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Hamiltonian 



k i=l 



The Fourier transformations are defined by 



air. 



_ | ifc(n-7/l 



-7/12-M/6) 



-ifc(n+7/12-m/6) 



(A3) 



where Z = 1,2,3 and m = 4, 5, 6. The Hamiltonian, up 
to the order of O(S), takes the form of 



H = E a + J2(uiai k a ik + -na ik af +lk + j-af k a i+lk ), 

k i=l 

where a 7k = a lk , E a = -2NS 2 (J AF - 2J F ), 

^1, 3,4,6 =S(JAF - Jf) + A*, ^2,5 =~2SJf + Mj 7l,2 = 

SJ F e~ % e, 73 = SJAFe~ l e, 74,5 = SJ F e l e, and 7 6 = 
SJafg 1 ® ■ Via the Bogoliubov transformation 

ajfe = «,i(fe)ai)t + u i2 (k)a 2k + u i3 {k)a 3k 

+uu(k)af k + u l5 (k)a+ k + u i6 (k)a£ k , 



jk 



= uji(k)aik + u ]2 (k)a 2k + u j3 {k)a 3k 

+u 3 ±{k)a+ k + u j5 (k)a+ k + u j6 (k)a+ k , (A4) 



where £ = 1,2,3 and j — 4, 5, 6, the Hamiltonian can be 
diagonalized as 



H = E g + ^ X a iEik a fk a ik, 
k i=l 

where cr^ = + for i = 1,2,3 and — for i = 4,5,6, 

and E g =^ + DE«f(E M a ) + X>*(E(Kf)]- 

fc i=l j— 4 2= 4 j — 1 

The coefficients of the Bogoliubov transformation can be 
found through equations of motion ihhi k = [an-,H]: 



determined through the constraint l|A2 

6 

6NS = EE n ^ 



k j=i 
3 



nik 



Tlmk 



i=l i=4 
3 6 

^ \u ml \ 2 (l + n ik ) + 2J I (A5) 



where Z = 1,2,3, to = 4,5,6, and Uj/s = {a1j k aj k )T- 
Defining hi k = (af k ai k ) F , the spin- wave distribution 
function, we have = [e°' i ' E "'' T — 1] _1 with <jj = + 
for i = 1, 2, 3 and — for i = 4, 5, 6. 

The internal energy and the magnetic susceptibility 
can be expressed as [l7| 



E = E, 



X 



y y &iEi k n ik , 

k i=l 
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For a given fc and /it, the eigenvalues E ik and eigenvectors 
{u\i, ...,v,6i) can be numerically calculated by the driver 
ZGEEVJ of the LAP AC K, which is available on the 
website [24^. 

At finite temperature, the Lagrange multiplier n(T) is 



where <Zj = + for i = 1, 2, 3 and — for i = 4, 5, 6. % 
and C within the LMSW framework can be thus calcu- 
lated numerically, as shown in Figs llOl 

For the PIMSW method, the HP transformation is 
treated, up to the order of 0(S°), such that 

^tn^2n — 2S , ai n a s J n — -a ln ai„ai n a 2 „ — -ai n a 2n a 2n a,2 n . 

The interactions may be handled in the perturbational 
way: 

a ii a in a i« a 2ri — 2(a ln ai„)odi n a2n + 2(ai„aJ n )oa^ u Oi n , 
a\ n a^ n a 2n a 2n ~ 2(ai„a 2n )oa 2n a2n + 2(a 2n a2„)oai«a2n- 

Define (a+ n a ln ) = C Uj (a 2 f n a 2 „) ee C 22 and 
(ai n a 2 h „)o = C\ 2 . We obtain 

G11 = jf E(4k^k)o - ^ E(l«i 4 | 2 + K 5 | 2 + \u 16 \ 2 ), 

k k 

c 22 = ^ E( a 2fc a 2/c) = j? E(l«2 4 | 2 + |m 25 | 2 + K 6 | 2 ), 

k k 

Cl2 = jj E( a lfe a 2fe)o = ^ E ( W H U 21 + "12^22 + ""13^23)^ 

fe fe 

and 

ai" n Oi n ffli n £^, = ^ 2(Ciiai fc a+ fe e -1 ^ + C 12 af k a lk ) 

n fe 
6 

= J2( 2 die' l!l uuu 2i + 2C 12 \ Uli \ 2 )af k a lk , 

k i=l 

E a l« a 2r l a 2 f n a 2n = X/ 2 (Cl2 a 2fe a 2fe + C 22 a lk a\ k e~ % ~Z) 
n fe 
6 

= EE( 2C ' 22e " <lu « u 2i + 2C 12 |'U 2i | 2 )a+a ife , (A6) 
fe i=l 
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where (...) in (|A6|) is the up-to-0(S , °) correction to Eik 
from the term (S^S^n), and the other terms have sim- 
ilar contributions to The only difference between 
LMSW and PIMSW is E lk , say, the former is treated 
up to 0(5 1 ), while the latter up to 0(S°). Accordingly, 
Eik, X an d C within the PIMSW scheme can thus be 
calculated numerically, as shown in Figs llUI 
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